A New Approach to Noninvasive-Prolonged Fatigue Identification Based on Surface EMG Time-Frequency and Wavelet Features

In sports, fatigue management is vital as adequate rest builds strength and enhances performance, whereas inadequate rest exposes the body to prolonged fatigue (PF) or also known as overtraining. This paper presents PF identification and classification based on surface electromyography (EMG) signals. An experiment was performed on twenty participants to investigate the behaviour of surface EMG during the inception of PF. PF symptoms were induced in accord with a five-day Bruce Protocol treadmill test on four lower extremity muscles: the biceps femoris (BF), rectus femoris (RF), vastus medialis (VM), and vastus lateralis (VL). The results demonstrate that the experiment successfully induces soreness, unexplained lethargy, and performance decrement and also indicate that the progression of PF can be observed based on changes in frequency features (ΔFmed and ΔFmean) and time features (ΔRMS and ΔMAV) of surface EMG. This study also demonstrates the ability of wavelet index features in PF identification. Using a naïve Bayes (NB) classifier exhibits the highest accuracy based on time and frequency features with 98% in distinguishing PF on RF, 94% on BF, 9% on VL, and 97% on VM. Thus, this study has positively indicated that surface EMG can be used in identifying the inception of PF. The implication of the findings is significant in sports to prevent a greater risk of PF.


Introduction
Surface electromyography (sEMG) is an electrical feld of human skeletal musculature [1]. It is acquired by placing electrodes on the skin surface near the human muscle. Te frequency and amplitude of the signals represent the behaviour and condition of the muscle's motor unit, conduction velocity, and ionic alteration of the muscle. Fatigue can be determined by the changes in its frequency content and amplitude either during an activity by analyzing every interval time length [2], at the beginning and ending of the activity [3,4], or before and after the activity [5][6][7][8].
In fatigue detection, frequency shifting represents the changes in muscle fbre conduction velocities and subsequent changes in the duration of the motor unit action potential waveform and fuctuations of muscle force and muscle fbre types as well as their decomposition [8,9]. Most of the opinions agree that fatigue can be identifed when its frequency shifts to a lower value to indicate that the muscle conduction velocities are slowing down [10,11]. Other than frequency, fatigue can be detected through the amplitude of sEMG signals. Te changes in the sEMG amplitude depend on the number of active motor units [12], discharge or fring rates, and the shape and propagation velocity of the intracellular action potential [10]. Te amplitude of sEMG tends to increase during submaximal voluntary contraction (during motor unit recruitment) and decrease during maximal voluntary contraction [10,[13][14][15][16].
Other than time and frequency features, a new timefrequency feature representation to track fatigue was introduced and is known as the wavelet index (WI) method [17]. WI was introduced since it is more suitable to deal with nonstationary signals such as sEMG [17]. Tere are fve WI features introduced by Malanda and Izquierdo, including the wavelet index ratio between moment −1 at scale 5 and moment 5 at scale 1 (WIRM1551), wavelet index of the ratio between moment −1 at the maximum energy scale and moment 5 at scale 1 (WIRM1M51), wavelet index of the ratio between moment −1 at scale 5 and moment 2 at scale 2 (WIRM1522), wavelet index of the ratio of energies at scale 5 and 1 (WIRE51), and wavelet index ratio between square waveform lengths at diferent scales (WIRW51). Trough WI, the distribution shifting of sEMG energy can be assessed based on its scale and frequency band of decomposition.
In normal conditions, fatigue usually disappears by itself after a while. Recovering from fatigue indicates that biochemical reactions during sports activity are able to return to a normal level [18]. Under normal fatigue (NF) conditions, most opinions agree that the degree of fatigue begins with an increment in amplitude, followed by unchanged and decreased trends, as well as accompanied by a decrement in frequency centers. WI features tend to increase, indicating the distribution of energy shifting to a lower value under NF conditions. However, high-intensity training activity will commonly lead to more biochemical reactions such as releasing of stress hormones (cortisol, epinephrine, and prolactin) [19,20], glycogen depletion [21], and the existence of lactate [22]. Fatigue due to intense training will require a longer recovery period than normal physical activity. It is crucial for improvement and recuperation [23]. During the period, it will enable hormones to return to a normal level [18] and allow physiological adaptation to a cardiovascular and muscular system to provide a higher level of performance [24]. If the training load is imbalanced with an inadequate recovery period, fatigue can be continuous and accumulated. Tis situation leads to prolonged fatigue (PF). Under this condition, more biochemical or maladaptive hormonal responses may occur [24,25]. Te alteration in biochemicals, which leads to PF, can be signifed by reduced performance, lethargy, soreness, insomnia, psychological disturbance, restlessness, hypertension, and increased incidence of injury [21,26]. It commonly requires several days to a week to recover from [23,27]. Tis condition needs to be treated accordingly to avoid a more severe condition, known as chronic fatigue syndrome. A report shows that about 20-60% of athletes, 60% of elite runners, and 33% of nonelite runners experienced chronic fatigue syndrome at least once in their career life [23,26,28].
In current practice, PF signs can be assessed invasively or noninvasively. Blood tests are invasive and used to investigate biochemical concentrations associated with PF such as lactate, glycogen depletion, creatine kinase, and iron levels [22,28]. Meanwhile, muscle biopsies are utilized to evaluate the condition of the injured muscle and ionic concentration in the muscle layer. During the collection of muscle tissue, numbing medicine is required. Although both blood tests and biopsies are reliable and accurate, they cause discomfort and are not suitable for frequent measurement. Furthermore, these methods are time-consuming, need to be analyzed in a laboratory environment, and require full supervision from an expert [29].
Due to the limitation of invasive methods, PF can also be traced through noninvasive diagnostic tools because the alteration in biochemicals can be observed physically. For example, glycogen depletion and lactate accumulation are commonly associated with a decline in performance, the oxidative stress increment leads to muscle pain, and cytokine leads to unexplained lethargy, decreased appetite, depression, and sleep disturbance [26,28]. Te commonly used noninvasive tools are interviews, athlete-coach monitoring approaches [18], questionnaires [26,28], training logs [30], and perceived exertion ratings. Te current practice requires more than one diagnostic tool to comprehensively screen of PF signs. Diagnostic tools such as interviews, training logs, and questionnaires often require close supervision by the practitioner and personal coach. Nevertheless, using many tools for the PF identifcation process is inefective, particularly, in monitoring a large group of athletes because these tools are time-consuming and have many procedures. Even so, many agree that PF condition prevention is the best solution [26]. Te reason is that the treatment of PF is timeconsuming and cost-inefective, depending on the degree of PF. Furthermore, PF signs endured are too risky for athletes.
Later, fndings reveal that the center frequency shifting of sEMG to the upper value was attributed to the alteration of ionic concentrations such as lactate and glycogen and the existence of soreness following high-intensity exercise [6,31]. Tis fnding is opposed to the earlier fndings that state a decrement in the center frequency of sEMG following short duration and light exercise refers to fatigue conditions. Tis situation demonstrates that duration, the intensity of exercise, biochemical reactions, and the existence of PF signs may afect the sEMG signal behaviour. Tis situation also demonstrates the potential of sEMG as a new tool which is noninvasive, comfortable, fast, easy to use, and quantifable to detect signs of PF. Te detection at the earliest stage helps prevent a more serious state of PF.
Terefore, this paper aims to investigate the ability of sEMG signals to identify the inception of PF in four muscles with diferent percentages of muscle activation. Tis paper also investigates the ability of wavelet index features in PF identifcation. Te performance of the surface EMG features was evaluated by the naïve Bayes classifcation accuracy in predicting the PF condition.

Study Protocol.
Twenty participants (age ± standard deviation (SD): 24 ± 3 years old; body mass index: 22.7 ± 2 kg/m 2 ) were recruited for this study. Participants were screened using a Physical Activeness Questionnaire (PAR-Q and You) (Supplementary Appendix S1) to rule out any pre-existing health contraindications and risk factors for exercise. Te exclusion criteria were participants with diabetes, high blood pressure, heart disease, any chronic disease, joint or bone problems, and taking any medication to control blood pressure and blood sugar. Te approval, to conduct the experiment procedure, was obtained from the Ethical Committee, Universiti Putra Malaysia (UPM/ TNCPI/RMC/1.4.18.1(JKEUPM)/F2).
Te participants were given a written and verbal explanation, including the potential risks and discomfort that they might experience. Te participants signed written informed consent before the experiment began. As a precaution, the participants were also protected by insurance (policy number: P809067176) during the whole experiment period.

2.2.
Procedure. Te experimental design was divided into two phases: Phase I was meant for familiarization and Phase II was for intensive training. Each participant had to take part in both phases. Phase I enabled the participants to familiarize themselves with the equipment and procedures, while Phase II was designed to induce PF signs. Between Phase I and Phase II, the participants were requested to rest and refrain from exercising or doing any heavy physical activity. Phase I was carried out on alternate days to avoid the emergence of PF, and Phase II was carried out on fve consecutive days (see Supplementary  Table S1. Schedule of Experiment).
Te participants were instructed to refrain from any heavy physical exercise, alcohol, and cafeine consumption 24 hours before the running test in both phases. Tey also required taking meals two hours before the assessment to avoid lack of energy and dehydration. Te experiment was conducted in accordance with the Bruce Protocol treadmill test (see Supplementary Table S2). In the protocol, the inclination and speed of the treadmill were increased every three minutes. Te total duration of the protocol was 21 minutes. Te participants were required to run for fve consecutive days and requested to improve their performance daily. As individual fatigue response is highly variable, no specifc distance and time duration were fxed [32].

Data Collection.
In this study, training logs (see Supplementary Appendix S2) were used to record measurements before, during, and after the running activity. Te measurements were used to monitor the daily performance and identify the emergence of PF conditions during Phase II of the experiment. Te fowchart of the experiment procedure and measurements is shown in Supplementary Figure S1, and the equipment utilized throughout the experiment was the COSMED T170 treadmill, Polar chest strap heart rate monitor, Watsons blood pressure monitor, and custom-made surface EMG data collection tool (see Supplementary Figures S2(a) and S2(b) for the schematic circuit of surface EMG systems).
Te recorded measurements during Phase II were as follows: (a) Percentage of the maximal heart rate Percentage of the maximal heart rate (%HR max ) is recorded to indicate running eforts performed by the participant. %HR max is determined as %HR max � HR max (running) X 100 (220 − Age) . (1) (b) Percentage of endurance time Endurance time of running on the treadmill is calculated [33] based on the following equation: %T endurance � T recorded X 100 21minutes .
(2) (c) Prolonged fatigue sign identifcation Te participants were also requested to fll in a 24-hour training distress questionnaire (see Supplementary Appendix S3) daily [28]. Te questionnaire was used to identify sleeping and psychological disturbance and muscle soreness. During the experiment, the participants were also interviewed before the running activity. Te interview session was conducted to identify whether the participants experienced lethargy. After running, the participant requested to scale the running activity experiment to indicate the difculty of the experiment. Te emergence of PF signs was monitored based on noninvasive diagnostic tools, as summarized in Table 1.
Te prolonged fatigue diagnosis was important as surface EMG signals were then grouped and classifed based on PF signs experienced by the participants. Due to ethical reasons and potential risks endured by the participants, only symptoms developed within fve days of the experiment were monitored. Te earliest PF signs that appeared during the training were sufcient to indicate the emergence of PF. Te participants were also reminded about two symptoms of fatigue. Te symptoms were observed based on two conditions as follows: (a) Fatigue symptom 1 (monitored before the running activity) Te participant was not allowed to run and was terminated from the experiment if any of the following fatigue symptoms were observed before the assessment: (i) Heart rate >100 bpm (ii) Blood pressure >140/90 (iii) Showing performance decrements in the previous experiment (iv) Psychology scores in the 24-hour training distress questionnaire >14 for at least three days (v) Collapsing in the previous experiment (b) Fatigue symptom 2 (monitored during running activity) Te participant must stop running if the following symptoms are observed while running: (i) Lack of energy (ii) Feel dizzy (iii) Blurred vision 2.4. Surface Electromyography. SEMG signals were collected from the biceps femoris (BF), rectus femoris (RF), vastus lateralis (VL), and vastus medialis (VM). Tese muscles were selected based on the activation muscles during running and sufer a high rate of injury in sports involving running [35,36]. Running at 10°grade inclination activates 79 ± 7% of BF, 76 ± 14% of vastus, and 44 ± 20% of RF, and the activation is elevated as the inclination increases [36,37]. Te sEMG signals were collected using the custom-built sEMG acquisition system, as shown in Figure 1.
Te sEMG system was designed based on an AD620 instrumentation amplifer system. AD620 was selected as it provides a 130 dB common-mode rejection ratio (CMRR), low power consumption, that is, 1.3 mA, and comprises a low input voltage noise of 9 nV/√Hz at 1 kHz and 0.28 μV p-p in the 0.1 Hz-10 Hz band. Te full system ofers a signalto-noise ratio (SNR) of 25 dB and gains an amplifer at 248. Te 50 Hz notch flter is designed according to the following equation: where R � 68 k Ω and C � 47 nF. Te schematic diagram of the sEMG data acquisition board is depicted in Supplementary Figures S2(a) and S2(b). Te analog signals of surface EMG were digitized into 12 bits by using National Instrument Data Acquisition (NI-DAQ) 6008 with frequency sampling, F s at 1k Hz. 1k Hz was selected to avoid aliasing as suggested by De Luca. Ten, the collected data were fltered using a digital fnite impulse response (FIR), a high-pass flter (HPF), 301 taps, and a cutof at 20 Hz. Tis HPF is essential for removing baseline wander during data acquisition. F s and HPF specifcations were set using data logger software, LabVIEW.
Ag/Ag Cl electrodes from Kendall MediTrace 200 were used to acquire the signals. Bipolar electrodes with 20 mm inner distance were attached to the involved muscles, and one reference electrode was placed at the knee of the participants. Te electrodes were positioned at BF, RF, VL, and VM based on the Surface EMG for the Noninvasive Assessment of Muscles (SENIAM) standard [38]. Te RF muscle was determined by 50% distance between the patella upper borders and the anterior iliac spine (AIS), VL was at 25% distance from Gerdy prominence to AIS, and VM was at 25% distance from the joint space to AIS. After measuring and marking the muscle, palpation of the involved muscle was carried out to ensure that the electrodes were placed correctly. During palpation, the participants were asked to fex and extend knee movements to activate the muscles, as shown in Supplementary Table S3 [38,39].
In data collection, the participants were asked to move their legs to activate the observed muscles. Only one leg was involved in data collection, and it was observed that the participants were comfortable using the left leg in the study. Te investigation on one leg was enough in this study to observe PF conditions based on surface EMG. RF, VL, and VL were activated when the hip was fexed, while the knee was extended to 180. As shown in Figure 2(a), the participants were asked to sit on a chair and were requested to move their legs from Point A to Point B to activate the quadriceps muscle group. Tey were asked to stay at each point for ten seconds and then repeat the movement three times. It was discovered that the RF, VL, and VM muscles contracted when the leg was at Point B and were at rest when the leg was at point A, as suggested by Konrad. Other than RF, VL, and VM, surface EMG signals were also collected from the BF muscle. Te location of the electrodes for BF was at 50% distance from the lateral epicondyle of the tibia to the ischial tuberosity. To collect surface EMG signals from BF, the participants were asked to stand and move their legs from Point D to Point E, as shown in Figure 2(b). Before that, the participants needed to place one (1) of their legs one foot (1 ft) away from Point C, which was Point D. BF contracted when the hip was extended, while the knee fexed [34]. When the body gesture was about 30 forward, the knee fexed until the leg was lifted about 15 cm to Point E. Tis distance was chosen as it provides the maximum activation of BF during movement [40]. Te participants were requested to move their legs from Point D to Point E three times at (10) seconds intervals. Te sEMG signals were collected during before and after running activities. Figure 3 shows the example of sEMG signals when the knee is fexed and extended during the position and movement in Figures 2(a) and 2(b).

Feature Extraction.
A nonoverlapping windowing technique was employed with samples n � 5000, as shown in Figure 4. Te moment at which the muscles started to contract and relax was ignored because the dynamic movement during data collection might result in false information [41]. Te number of n was selected because the authors of [2] have demonstrated that segmentation length is suitable for muscle fatigue identifcation. Te features were extracted at each contraction and averaged.
Specifc features were extracted based on the frequency, time, and wavelet index (WI). Te spectral content of surface EMG was determined according to the Fourier transform, and the frequency parameter was quantifed based on its median (F med ) in (4) and mean (F mean ) in (5): Fmean where f j � frequency of the spectrum at frequency bin j, P j � EMG power spectrum at frequency bin j, and M � length of the frequency bin. Te time features can be quantifed based on the mean absolute value (MAV) in (6) and the root mean square (RMS) [42] in (7) [42]: where x = signals and n = number of samples [43]. Tis study also investigates the ability of WI features in PF [43] identifcation since it was never tested in determining fatigue under high-intensity conditions. WI features were used to evaluate the distribution shifting of sEMG energy based on its scale and frequency band decomposition. WI was calculated based on the discrete wavelet transform (DWT) which was decomposed into fve levels by using symlet 5 (sym5) and Daubechies (db5) as the mother wavelet [17]. Te decomposition process consisted of a series of flter banks, where at every i level of decomposition, the signal was fltered into half of the frequency band [44]. Te low-pass flter produced an approximation coefcient, while the high-pass flter produced a detail coefcient (D i ) (scales). Figure 5 shows the fve levels of sEMG decomposition details and the power spectra of decomposition details at scales 1-5 that are determined based on the Fourier transform.
Te wavelet index ratios between moments at diferent scales were then determined based on the power spectrum of wavelet details, D i. .
Te fve WI features were tested as follows: (a) Te WI ratio is between moment −1 at scale 5 and moment 5 at scale 1 (WIRM1551).
where ψ(t) used was sym5, f1 = 10 Hz and f2 = 500 Hz, and D 5 (f ) and D 1 (f ) are the power spectra of the fve and frst scales of decomposition details [14]. (b) Te WI ratio is between moment −1 at the maximum energy scale and moment 5 at scale 1 (WIRM1M51).
During the extraction of WI features, frequency sampling F s � 1k Hz [44] and n � 1024 were used. Te WI features were then log-transformed to follow the normal distribution.
Fatigue identifcation always refers to the increment or decrement in the features before and after the activity [45]. Te changes and shift of the features (F) in this study were quantifed by Δ F � Fpost -Fpre.
Te positive value of ΔF indicates a feature increment for postexercise, whereas the negative value indicates feature decrements. were distinguished based on PF signs explained in Table 1. While the features of the participants who did not experience PF conditions were grouped into NF, the features of the participants who experienced PF conditions were grouped into PF. A t-test was conducted, and a signifcant value was set at P < 0.05.

Daily Plot of Surface EMG Behaviour.
Te daily plot of sEMG behaviour for NF and PF conditions was performed to investigate the progression of fatigue in diferent muscles with diferent activation percentages. It was plotted based on ΔF med and ΔRMS since these two features commonly represent time and frequency information on surface EMG in   Journal of Healthcare Engineering fatigue identifcation. A daily plot was also conducted based on fve ΔWI features. To identify the changing behaviour of the features during the emergence of PF, the features were normalized by plotting them from a day before the emergence of PF and the frst three days under PF conditions. Te reasons were the individual's responses to PF signs that varied, and the normalization will help understand the trend line of sEMG features during the intensive training period specifcally under PF conditions. From the feature selection, dimensionality reduction was employed to reduce the complexity and computation time of the classifcation algorithm, increase accuracy, and decrease overftting problems [46]. Data reduction in this work was carried out based on linear discriminant analysis (LDA). Tis method maximizes the intercluster distance between classes and minimizes the intracluster distance within classes in the transformation of reduced features. In LDA, the original dimensional feature space was transformed into a lower dimensional feature space, without losing any important information [46].
In the classifcation stage, the naïve Bayes (NB) technique was applied to discriminate NF and PF classes. Tis method was selected as it was previously utilized in experiments studying fatigue classifcation [36,37]. NB is one of the established statistical pattern recognition methods [46]. NB classifer functions are based on the probability distribution of the feature vector, x. x belongs to class ω m which is computed from probability distribution conditioned on the class ω m, P(x|ω m ), by assuming class-conditional independence of the features: where d is a dimension of the feature instance x. Equation (13) requires that the k-th features of the instance, which is x (k) , are independent of all other features, given the class information.
Te probability of the x class itself is characterized by Te classifcation rule was computed from the discriminant function g m (x) to represent posterior probabilities as It was represented for each m-class. Meanwhile, the x class is determined by the largest g m (x) computation.
In this work, k-fold cross-validation (CV) was adopted for training the classifer. Te performance of the classifcation was evaluated for the accuracy, specifcity, precision, and average CV error (CVErr). Table 2 shows a daily % HRmax record. It indicates that about 18 participants ran at their maximal efort by showing %HRmax >80%, based on the Edwards Intensity Zone 1992. Running at this rate caused the participants to experience heavy breathing and muscular fatigue. It proves that the Bruce Protocol treadmill test provides the high training intensity required in this experiment. High-intensity exercise is essential for inducing faster PF signs. Physiological fatigue responses under PF conditions are tabulated in Table 3. It shows that the frst PF sign developed was muscle soreness, which was on day 2 (D 2 ) of the assessment. Tis situation was found to be similar to other studies in [6,47], whereby soreness developed as early as 24 hours after strenuous exercise. Table 3 also indicates that PF signs accumulated with performance decrement starting at day 4 (D4) of intensive training. Te results agree with [16] as the untreated PF condition develops more PF signs. Apart from that, the results suggest that only three PF signs appeared within fve days of intensive training including soreness and performance decrement. Moreover, the results reveal that these are the earliest signs of PF developed in the study. Te result in Table 3 further shows that none of the participants experienced psychological and sleeping disturbance, restlessness, and hypertension following intensive training. Hence, the classifcation of collected surface EMG signal features was based on physiological responses identifed in Table 3. Te term PF condition afterward refers to muscle soreness, performance decrement, and lethargy.

Surface Electromyography.
Te daily plot of sEMG feature behaviour is displayed in Figure 6, while the bar plot represents standard deviation, "o," and "x," symbols represent the mean value of features in NF and PF conditions, respectively. Te features under NF conditions were plotted from day 1 to day 5 (D 1 -D 5 ) of the assessment, whereas the plots under PF conditions were normalized from the day before the emergence of PF (D NF ) to the frst three days of the occurrence of PF signs (D 1 -D 3 ).

Frequency Feature.
Teoretically, the frequency information changes in sEMG describe the behaviour of conduction velocities inside the muscle and subsequent changes in the duration of the motor unit action potential waveform and fuctuation of muscle force and muscle fbre types as well as their decomposition [8,9]. Te frequency spectrum shift information is represented by its mean (F mean ) and median (F med ) in assessing muscle fatigue [42]. Figure 6 shows that ΔF med resulted in a negative value for BF, RF, VL, and VM under NF conditions. Te negative values of ΔF med demonstrate that F med was decreasing postrunning activities. Te decrement in F med was like the most dominant opinion where frequency tends to shift to a lower value to characterize fatigue. Te decrease in the centre of frequency as a result of reduced muscle conduction velocity and a change in the frequency spectrum was brought on by the absence of high threshold motor unit recruitment. However, ΔF med shows positive values for BF, VL, and VM on day 4 (D 4 ) and day 5 (D 5 ) under the NF condition. Te positive values of ΔF med indicate that the median frequency spectrum was shifted upwards. An increase in F med was also identifed on the day before PF signs appeared (D NF ) for the BF, VL, and VM muscles, and this behaviour was sustained throughout the PF condition.
Te plot in Figure 6 also indicates that the positive values of ΔF med only occurred in RF during PF conditions. Statistical analysis reveals that an increment in ΔF med under PF is signifcant at P < 0.05 for BF, RF, VL, and VM, as tabulated in Table 4. An increase in F med of sEMG during fatigue was rarely reported. Te increasing center of frequency was once reported by [48] during the frst 30 minutes of recovery from dynamic exercise at a load of 80% of the VO 2 max. Te increasing center of the frequency was reported following the elevation of temperature and lactate after high-intensity dynamic exercise [48]. Te relationship between the skin and muscle temperature and the increasing center of the sEMG frequency spectrum was later confrmed in [49]. Te positive linear relationship between the temperature and median frequency might be due to an increase in the muscle conduction velocity to increase the power output [49,50]. Te relationship between the frequency of sEMG and temperature was also discussed in [51].
In [51], the authors demonstrated that there is a less efect of temperature on muscle strength and frequency of sEMG but related other possibilities that afected the frequency features such as diferent recruitment properties of the motor units and the percent of fats and slow twitch motor units under electrodes. Heavy dynamic exercise might also contribute to the substitution of muscle groups following an efect on the alpha motor neuron pool through refex inhibition that alters recruitment properties. Te efect of the neural drive on the muscle and its motor unit action potential was also identifed as one of the factors that afect sEMG components [52].
Te neural drive for the muscle factor might be related to fatigue induced in the peripheral and central systems. Fatigue in the central system occurs when neurochemical in the brain is altered and stress hormones are secreted. When this happens, it will modify the peripheral information in the contracting muscles and afect the characteristic of sEMG [53][54][55][56].
Te increment in frequency information features ΔF med under PF conditions also might be due to fatigue at the peripheral system. Fatigue at the peripheral system arises from the muscle itself when there is impairment of the peripheral mechanism due to high-intensity exercise as demonstrated by participants in this experiment [45,49,50]. High-intensity exercise reduces blood fow due to intense muscle contraction which causes the inadequacy of oxygen supply to the muscle. Tis situation is also known as an anaerobic condition [54]. Te inability to get enough oxygen triggers a biochemical reaction in allowing muscle contraction [57,58]. An inadequate recovery period causes the inability of ionic alteration during high-intensity exercise to return to its normal level and continue to accumulate. Tis situation is signifed by the emerging PF signs such as soreness and performance decrement. Te ionic changes most probably involve glycogen breakdown and the presence of lactate concentration. It is supported by the recorded %HRmax during the running activity, of which 80% and above commonly involves anaerobic contraction. In anaerobic contraction, glycogen and lactate concentration play important roles in ensuring muscle continuous contraction [57,58]. Furthermore, the alteration in glycogen stores normally leads to soreness and performance decrement due to inadequate fuel for workload [22,31], and the release of lactate contributes to fatigue and muscle pain, as experienced by the participants in this study. Tis situation is supported in [31,48] that also demonstrated that the alteration of both concentrations led the frequency of sEMG to shift to the upper value. Figure 6 also demonstrates that there were diferent increment trends in ΔF med among the investigated muscles. Te trends that happened might be related to muscle activation during running activity. Running at a higher slope such as in the Bruce Protocol treadmill test requires more muscle activation from BF, VL, and VM than from RF, as demonstrated in [36,37]. A previous study shows that more muscle activation leads to faster progression of fatigue [59]. Tis study has demonstrated that changes in ΔF med happen faster in more activated muscles than in less activated ones. Te fast changes in the frequency made the observation of PF signs through more activated muscles rather difcult. Te reason was the frequency feature increased even without the emergence of PF signs. Indirectly, an increase in the median frequency might be due to an increase in the muscle temperature and muscle conduction velocity to increase the power output, substitution of the muscle group and recruitment properties, and alteration of ionic concentration underlying the muscle that progressed faster in muscle activation during running. Tis study has also proved that PF conditions could be easily observed from less activated muscles such as RF because the increment in frequency only occurred under PF conditions. It indicates that PF can be easily identifed when frequency from less activated muscles starts to increase.

Time Features.
Muscle activity can be observed through its amplitude during the contraction in time-domain representations. In fatigue identifcation, changes in its amplitude signify the degree of fatigue experienced by the subjects. As exhibited in Figure 6, ΔRMS of BF, RF, and VM under NF conditions increased on D 1 and decreased on the following days. Tis progression is similar to dominant opinions that with an increment in amplitude, the decrement in behaviour in characterizing the degree of fatigue soon follows [10].
However, the ΔRMS increased again, as shown in D5, in the BF and VM muscles ( Figure 6). Teoretically, in normal conditions, when the load increases, the amplitude tends to have a larger decrement. In this study, the load refers to the endurance time, for which the participants were asked to improve their performance daily. Based on the plot of ΔF med on similar days D 5 on BF and VM, the median frequency shows an increment. In the previous section, the increments in ΔF med were related to an increase in temperature. However, the fndings in [48,49] have shown that the increment in the muscle and skin temperature will reduce the amplitude of sEMG signals. Te increment in ΔRMS indicated by the BF and VL muscles in D 5 might be due to the release of free-resting calcium which resulted in force potentiation and led to the increment in EMG, as demonstrated by [60].
Te increment in ΔRMS especially under PF might also be due to the changes in ionic concentration. Te changes in the ionic concentration were observed through frequency feature behaviour in the precious section. Te Journal of Healthcare Engineering fndings in [61] reported that there was a curvilinear positive relationship between lactate concentration (after reaching a certain lactate threshold) and the amplitude of sEMG. Apart from that, Figure 6 discloses that ΔRMS started to decrease again on D 2 of PF, specifcally in BF and VL. Te decrement in ΔRMS under PF conditions was discovered in [31,62]. Both studies have proved that amplitude decreases during the emergence of soreness. Nevertheless, another study reveals that the amplitude increases under similar conditions [63]. Terefore, it is reliable to state that the amplitude increases or decreases under PF conditions. Te increment and decrement in amplitude under PF also show the degree of fatigue experienced by the muscle. Tis is attributed to the decrement in amplitude under PF which occurred in highly activated muscles like BF, VL, and VM. High activation led to the fast progression of fatigue. It began when the frequency features started to increase, followed by the amplitude which also increased. Ten, it continued with the decreased behaviour to show a certain degree of fatigue experience. Tis fnding was also supported by the progression of fatigue mapped on RF, which was less activated in the study. Figure 6 (Table 4). Te statistical test also indicates that the behaviour of ΔRMS under both conditions for BF, VL, and VM is not signifcant at P < 0.05 due to the fuctuation trend in the daily plot (see Figure 6).

Wavelet Indices.
Tis study also investigates the ability of WI features in PF identifcation. Te fve WI features were studied, as proposed by [43]. WI features tended to have similar behaviour and response to BF, RF, VL, and VM, as observed in Figure 6. Tey also tended to increase under NF conditions and decrease under PF conditions. Figure 6 for ΔWIRM1551, ΔWIRM1M51, and ΔWIRM1522 illustrates the transition of the increment and decrement in WI features under NF conditions. It is also important to note that the features were constantly decreased under PF conditions. Te increment (positive value) in WI features under NF in Figure 6 is similar to the results demonstrated in [17]. Te increment in the features specifes that the energy distribution shifted to a lower frequency band indicating similar behaviour of frequency, which tended to decrease to show fatigue conditions [10]. Figure 6 also demonstrates that the increment and decrement transitions occurred faster in high-activated muscles such as BF and VL. Tese situations can be observed under NF conditions on D 4 and D 5 for ΔWIRM1551, ΔWIRM1M51, and ΔWIRM1522. ΔWIRM1551, ΔWIRM1M51, and ΔWIRM1522 features also demonstrate that PF could be easily identifed in RF, as it occurred on the ΔF med and ΔRMS daily plot. Te decrement in ΔWI features was caused by energy distribution which slowly shifted to a higher frequency band, which caused the energy distribution at the lower frequency band of decomposition to decrease.
Te WIRE51 feature was quantifed in accord with its coefcient details D of decomposition. Te increment in ΔWIRE51 features showed a higher value of D at level 5 postexercise than preexercise. Figure 6 indicates that similar trends also appeared in another daily plot of WI features, of which ΔW IRE51 gradually decreased under NF plots. Furthermore, the value constantly decreased under PF conditions. Although D indicates the time representation of decomposition, ΔWIRE51 proved that the behaviour of the features did not rapidly fuctuate as demonstrated by ΔRMS behaviour. Te robustness and sensitivity of WI in dealing with nonstationary behaviour in sEMG were exhibited.
ΔWIRW51 was used to show accumulated changes in the waveform length ratio at D level 5 to D level 1. Trough waveform length behaviour, the duration, frequency, and amplitude of the surface EMG signals were efectively compressed [17]. Te increment in ΔWIRW51 in Figure 6 under NF suggests that the surface EMG waveform fuctuated faster postexercise than preexercise. Te features gradually decreased to indicate the fuctuation of the surface EMG waveform at D level 5 which was getting slower during postexercise. ΔWIRW51 persistently decreased under PF conditions. Hence, it signifes that, apart from the amplitude and energy distribution in the spectra, the waveform characteristic of surface EMG also changed due to PF.
Although BF, RF, VL, and VM demonstrate similar behaviour of WI features under NF and PF conditions, statistical results indicate that all fve ΔWI features are only signifcant at P < 0.05 for the RF muscle, as tabulated in Table 4. Table 5 indicates the classifcation results based on the NB method in identifying PF conditions. Tis result reveals the ability of sEMG features to distinguish between NF and PF based on the naïve Bayes (NB) classifcation method. Table 5 tabulates the lowest accuracy results from time features of BF, VL, and VM, due to fast fuctuation and overlapping plots of time feature values displayed in Figure 6. Tis condition makes predicting PF conditions through these features quite difcult. Te results revealed in Table 5 indicate that the frequency features had better classifcation accuracy than time features. Better accuracy was assisted by the signifcant statistical test results and daily plots to distinguish between NF and PF of frequency features. Te results in Table 5 reveal that the feature selection based on time and frequency ofers high-performance accuracy, specifcity, and precision in comparison with other feature selections. Tus, it can be concluded that both the time and frequency features of sEMG are signifcant for PF identifcation. In this study, the combination of time and frequency feature selections ofers accuracy at a rate of 94% on BF, 98% on RF, 95% on VL, and 98% on VL in distinguishing PF conditions. Te classifcation of performances in Table 5 proves the ability of WI features in PF detection. Te result shows that WI features produced good classifcation accuracy in BF (82%), RF (91%), and VL (80%) and less in VM (66%).

Conclusions
In conclusion, this study has demonstrated that the presence of PF can be identifed using the surface EMG signals. Te study also introduced a new quantitative noninvasive method to monitor the progression of fatigue, specifcally in the muscle of athletes. Tis monitoring method can provide information to athletes on their performance, and they can perform at their optimum energy. Tis noninvasive method is suitable to be applied in the sports feld for fatigue management and prevent chronic fatigue syndrome for athletes.

Data Availability
Te surface electromyography physiology data used to support the fndings of this study have not been made available because they involve the third-party right and participant privacy.

Conflicts of Interest
Te authors declare that there are no conficts of interest regarding the publication of this paper.

Acknowledgments
Tis work was supported by Universiti Putra Malaysia under IPS Putra Grant.

Supplementary Materials
Table S1: schedule of the experiment. Table S2: Bruce Protocol treadmill test. Table S3: prolonged fatigue sign identifcation. Table S4: quadriceps muscle movement based on the fexed and extended knee. Appendix S1: PAR-Q and You. Appendix S2: training log data collection form. Appendix S3: 24-hour history training distress questionnaire. Figure S1: fowchart of the experimental procedure and data collection. Figure S2: schematic circuit of the surface EMG data acquisition system. (Supplementary Materials)